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Abstract 

Several models of quantum open systems are known at present to violate, 
according to principles of the standard quantum theory of open systems, the 
second law of thermodynamics. Here, a new and rather trivial model of an- 
other type is suggested describing mechanism that violates, according to the 
same principles, the zeroth and the second laws of thermodynamics simulta- 
neously. Up to a technically minor modification, the model resembles some 
models already known, solved by standard means, and properly understood. 
Universal validity of two basic principles of thermodynamics in strictly quan- 
tum situations is thus simultaneously called in question. 
PACS numbers: 05.30.-d, 05.70.-a, 44.90.+C 
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I. INTRODUCTION 



Phenomenological arguments against general validity of standard statistical thermodynamics 
and call for inclusion of cooperative, selforganization and similar complicated phenomena even 
in absence of, e.g., external flows and far from equilibrium exist already for a long time. Since 
early nineties, one can find them especially in, e.g., theory of electron-transfer chemical reactions 
where, in connection with their phenomenological non-linear description, inclusion of such effects 
seems to be indispensable Necessity of their inclusion follows also from detailed analysis of 
what is known in molecular biology about how individual molecules (molecular machines) work in 
living organisms Q. Selforganization is usually believed to be a domain of nonlinear theories. In 
1996, the first hamiltonian linear quantum model was suggested that can lead to a selforganized 
state upon thermalization in a bath even when no external flows exist and, simultaneously, this 
state is energetically disadvantageous HQ . Recent analysis has revealed that the former non- linear 
phenomenological and the latter linear first-principle type of reasoning strive in the same direction 
and can be easily united The reader is referred to |6j-||] or also to || for previous models 
(extending also reasoning of ||,||]) where the cooperative and selforganizational tendencies in such 
models can be shown to provide a basis and unique possibility of violating even the second law of 
thermodynamics. 

Some of such models allow a mathematically rigorous treatment throughout all the calcula- 
tions ( H is perhaps the first one of them). However, for technical reasons, just one-step processes 
have been treated so far. Only very recently, first in 1998, the first rigorously solvable quantum 
model working cyclically as a perpetuum mobile of the second kind (i.e. converting heat from a 
single bath into, this time, a usable work without compensation) and violating thus, for the first 
time explicitly, the Thomson formulation of the second law of thermodynamics || was reported 
|lO|-|l^|. In 1999 and 2000, also other groups arrived, independently and for other situations, at 
the same conclusion challenging universal validity of the second law [13|-15]. From them, in par- 
ticular paper [|l4|] by Allahverdyan and Nieuwenhuizen inspired a public response |l6| . Completely 
different mechanism (connected with dynamically maintained steady-state pressure gradients in 
rarified gasses) potentially allowing violation of the second law was recently suggested by Sheehan 
- compare, in connection with previous paper [17], the discussion in |lg|,19[|. Another and even 
rather positively experimentally tested system was suggested by Sheehan already in 1994 [20-22]. 
Because of complicated nature of the problem as well as owing to 150 years of traditionally pre- 
sumed universal validity of thermodynamic principles, it is likely that irrespective of final result of 
the above counter-examples, the problem of potential violation of thermodynamic principles will 
remain topical for many years to come. A comment is only worth mentioning here that the above 
models of the present group always work from strong correlations (entanglement) among particles 
and/or competing and mutually interfering reaction (transfer) channels. Thus, like in l^,14j], the 
mechanisms discussed are appreciably different from those based on, e.g., the Feynman ratchet and 



pawl systems [53,24]. It is most likely that such ratchet-like systems (that would transform the 



thermal noise into one-directional linear or rotational motion) rather fail in experiment once they 
are devised to violate the second law [24], irrespective of previous opposite expectations but fully 
in accordance with the Feynman original analysis p3[ . 

The above models of the present group (usually described as isothermal Maxwell demon models 
because of dynamic opening and closing 'gates', i.e. reaction channels) are in fact so far sufficiently 
complicated. Moreover, their classical counterparts do not work. All this is why they can and 
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really do, at the first inspection, naturally induce fully comprehensible prejudice, mistrust, or 
misunderstandings. To this mistrust, also the fact contributes that the principles on which their 
activity relies remind of the original Maxwell demon |25| (opening and closing a gate after checking 
state or performance of previous steps and thus deciding about next elementary steps in, e.g., 
particle transfer). The idea of Maxwell demon was, however, often and in detail analyzed during 
the last 130 years ^f|^7]]. Result of the analysis was usually negative but it well applies to just 
classical models. In this connection, it is worth mentioning that all models P^]i~C|~14| really 
cease to work in the classical (high-temperature) limit. On the other hand, this analysis does 
not regard purely quantum models working on, e.g., principles known from nature, in particular 
from the contemporary molecular biology (interplay between particle transfer and accompanying 
topological reconstruction of the particle surroundings built in [j^— |^JTo|— 12. 

The situation seems to be even more serious as the quantum model of the isothermal Maxwell 
demon of (model of uni-directional isothermal particle transfer even against potential forces) 
allows a simple generalization to a greater set of sites available to a greater number of particles. This 
analysis then implies that in a stationary state (equilibrium one in the sense of thermodynamics), 
chemical potentials of one sort of particles in two subsystems interconnected by sophisticated 
(e.g. molecular) bridges could become even different pgj| . This questions universal validity of 
another basic principle of the statistical thermodynamics ||29[|. In order to make the situation 
simpler, we have here rebuilt the model of |3],|| to describe transfer of excitons (i.e. excitation 
energy) and analyzed principles of its work in connection with those of the original model. It 
appeared that the rebuilt model could be appreciably simplified to such an extent that it becomes 
fully independent. No possibility is seen to simplify it further and to preserve, simultaneously, 
the unusual phenomena investigated. The above elementary steps of 'checking performance of 
previous steps' and 'opening or closing the gate according to the result of the check', so typical 
of the Maxwell demon - like models, completely disappeared. What, on the other hand, remained 
is the existence of quantum interference of different reaction (transfer etc) channels. Except for 
a small (rather physical than technical) modification connected with presumed initial conditions 
and existence of two thermodynamic baths, the form of the model is fully standard. Without this 
modification, its solution is known, has been obtained many times and in many different ways, and 
is correspondingly believed to be well understood. Simplicity of the model and many times verified 
applicability of the really standard technical 'weaponry' applied to it is what may then, hopefully, 
change the so far reserved attitude of general public to the above provoking ideas questioning, on 
grounds of the quantum theory of open systems, universal validity of principles of the statistical 
thermodynamics. These questioned principles include now, in the light of the present results, not 
only the second but, as argued below, also the zeroth law of thermodynamics. 



II. MODEL 

System of a few levels interacting with a thermodynamic bath is a standard quantum problem. 
We shall use it also here. Specifically, here, we assume three levels, one ground and two excited 
ones, and refer to the excited levels as those with a Frenkel exciton placed either on site 1 or site 
2. These sites might be, e.g., two different molecules or local centres in, possibly, two different 
but adjacent solids representing two different electron subsystems. As it is easy to verify, it is for 
our problem here irrelevant whether we include or ignore the fourth level corresponding to both 
molecules (molecular systems, local centres) excited. For technical simplicity, we choose the latter 
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alternative. Corrections owing to the (here ignored) two-exciton states are, in, e.g., the excited 
state occupation probabilities and for small 1 <-> 2 exciton transfer rates, oc exp{— j3{e\ + £2)} where 
€j, j = 1,2 are the local exciton energies. Thus, they can be easily identified. 

The exciton residing possibly at sites 1 or 2 (if not lacking at all for some time owing to 
finite life-time effects admitted here) can in principle be transferred between the sites (subsystems) 
either coherently or incoherently. Here, we choose the first alternative, designating the hopping 
(resonance or transfer) integral as J. Finally, we complete the model by adding a bath interacting 
with the system. The tricky feature whose real sense will be seen only below is that we ascribe 
to each subsystem (designated as I with site 1 and II with site 2) its own thermodynamic bath 
represented by harmonic oscillators (phonons). Because we assume that initially, both the baths 
have a canonical distribution with possibly equal temperatures, this step is in such a case isomorphic 
to assuming that the exciton at site 1 or 2 interacts respectively with, e.g., just even or odd modes 
of a single bath. Technically, the next step is not relevant but physically, it is important to assume 
that these two baths are in a way separated as we shall discuss the energy (heat) flow from bath 
I to II and vice versa. Adding that the exciton-phonon coupling (with coupling constants g K and 
G K ) leading to the above exciton finite life-time effects (i.e. non-preserving the number of excitons) 
is assumed linear in the phonon creation (bj K ) and annihilation (bj K , j = 1 and 2) operators, one 
can directly write the Hamiltonian as 

H = Hj + Hjj + J(a\a 2 + aim), 
Hi = eia\ai + fno K b\ K bi K + -= g K fLU) K {ai + a\)(b lK + b\ K ), 

K ViV K 

Hn = e 2 a\a 2 + huj K bl K b 2K + — L= ^ G K fiuj K (a2 + 4)(6 2k + 4k) • i 1 ) 
k \N K 

The exciton creation (annihilation) operators are assumed to fulfil the Pauli relations 

{ai,a|} = {a 2 ,4} = 1, [ai, a\] = 0, 



{at, ai} = {a 2 , a 2 } = [01, a 2 ] = etc. (2) 

Here, {...,...} and [...,...] are the usual anti- and commutators. The phonon frequencies u K are, 
for simplicity, assumed the same for both the baths. Finally, €j, j = 1, 2 are the exciton energies 
while iV is the number of the phonon modes (finite before taking the baths thermodynamic limit) 
in each bath separately. Technically, the existence of two separated and uncorrelated baths, each 
of them interacting with just one exciton level, makes also the form of our matrices below simpler. 

We now want to proceed by writing down a closed set of equations for the exciton density matrix 
only, projecting off the information about baths. There are several well known ways leading finally 
to the same result. In order to be specific, we choose time-convolutionless Generalized master 
equations []30|-[3^| with the Argyres-Kelley projector [^J3^| . Practical application of the resulting 
equations then requires approximations (e.g. expansions) upon calculation of coefficients in the 
equations that can be avoided by application of scaling arguments. The method is well developed 
and standard now. Describing the formal apparatus, we must be, however, more specific as we are 
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now going to deviate from a standard weak-coupling approach. This is necessary here because of 
necessity to describe properly interplay between/among several competing processes. 

The weak coupling theory was made mathematically perfect by Davies who properly implied 



the scaling idea of van Hove [37.38]. It is funny to realize that, at least for systems of finite number 



of levels, the mathematical prescription how to calculate properly the weak-coupling dynamics 



as provided by the Davies theory is not unique. For that, compare Theorem 1.4 of [38] which 
establishes full uniqueness just in the full Van Hove limit [see (Q) below] but not for arbitrarily 
weak though finite coupling strengths as it corresponds to reality. Irrespective of it, because of its 
mathematically rigorous form, this theory (really rigorous in its region of validity) so influenced 
theoreticians that many of them now consider the weak coupling language as universal. This is, 
of course, unintentional negative consequence of the mathematical precision of the Davies theory. 
So, it is also (and even more than above) funny to observe that this theory fails to describe, e.g., 
the physical regime we are here interested in. This is owing to the overestimation, in the Davies 
(or, more generally, weak coupling) theory, of the role of in-phasing (that is, in our case, owing to, 
e.g., the coherent exciton transfer 1 «-> 2) as compared to the standard dephasing. [] The latter 
dephasing process is, in our case, simply due to the exciton (electron) coupling to the bath. Let 
us be even more specific: 

The standard weak-coupling theory is based on the notion of a small parameter, say g, of 
the coupling of the system to the bath. This means that g would be a joint small parameter of 
the last terms on the right hand sides of Hj and Hjj in (||) only. Then a new unit of time, say 
t = t/t' is chosen and we work, instead of the true physical time t, in terms of the dimensionless 
time if . Avoiding here technical details how to avoid Poincare cycles by performing first the 
thermodynamic limit of the bath [39], the result is that finally, the proper weak-coupling equations 



for the exciton density matrix are obtained (otherwise as below) by taking the combined limiting 
Van Hove procedure 

g — ► 0, t — ► +oo, g 2 T = const. (3) 



Omitting at this moment also less important mathematical details (see |3qj40[ ]), the statement is 
that the resulting equations then describe properly the time development of the system, in terms 
of the new time t' , to its canonical state corresponding to the (initial) temperature of the bath. 
Because of the limit g —> incorporated into @ and because the coherent (transfer or hopping) 
integrals responsible for the in-phasing are kept finite during this limit, the role of the dephasing 
processes in competing the in-phasing ones is fully suppressed. This is why the the corresponding 
asymptotic, i.e. canonical density matrix is then diagonal in the basis of the eigenstates of the 
Hamiltonian of the system H$ alone. In other words, this is why the relaxation goes to eigenstates 
o£H s . 

In order to incorporate the above competition between dephasing and in-phasing processes, we 
shall proceed in almost the same way except for one point: We take g as a small parameter of not 
only the system coupling to the bath. We also assume that J oc g 2 . The resulting limiting procedure 
is thus not that of the weak coupling but that of, rather, slow transfer processes. Concerning the 



This overestimation is not owing to any formal error, it is because of the very form of the 
presumed Van Hove scaling. In Nature, there is, e.g., no possibility to scale coupling constants as 
they are real constants, not variables. 
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relative strength of the 1 <-> 2 transfer and relaxation owing to the coupling to the bath, it might 
be, in such a scheme and for a general situation, still arbitrary. Here, we shall, however, assume 
that the latter coupling is rather intermediate or even strong as compared to the coherent 1 <-» 2 
transfer of the exciton in the sense that the dephasing is either comparable to, or even dominates 
over in-phasing. (Concerning the words '...dominates over..' remember, that in the situation of 
the weak coupling case, the in-phasing is, in the limiting sense of (|3|), infinitely stronger than the 
dephasing. Here, on the contrary, we admit in our case of the intermediate or strong coupling 
that the dephasing could be even dominating over the in-phasing but, figuratively speaking, their 
ratio remains always finite though perhaps arbitrarily large.) This causes remarkable differences in 
structure of the equations we aim at as well as in the physical conclusions. In particular, we then 
get that the relaxation does not go (in the sense of diagonalizing the asymptotic density matrix) 
to the eigenstates of the Hamiltonian of the system. To what state (i.e. to what exciton density 
matrix) the relaxation then goes we shall see later. 

Technically, the method of deriving the closed set of equations for the exciton density matrix 
proceeds in the following steps: 

• We introduce the density matrix of the 'system+bath(s)' complex in the 'interaction' picture 
as 

p(t) = exp{i£ • t}p s +B(t). (4) 



Here ps+B(t) is the density matrix of the system and the bath in the Schrodinger picture 
and the Liouvillean Cq ... = [Hq, . . .]/% where, however, Hq = J2'j=i[ e j a j a j + J2 K ^K^jjjj^. 
In other words, the hopping term J(a\a2 + a 2 ai) is now not included in Hq. This is unlike 
the scaling inherent to the weak-coupling case (0). 



We apply, e.g., the Fulihski and Kramarczyk identity 3C,31|, or its more famous form by 



Shibata, Hashitsume, Takahashi, and Shingu [32, 3^] 

jVp(t) = -iVC(t)[l + ij\xp^{-ij\l-r)C(T 2 )dT 2 }(l - V)C(n) 



•exp_^{z / C(t 2 ) dr 2 } dri) 1 
JtI 

[exp^{-i f\l - V)}C(t) dr}(l - V)p(0) + Pp(t)]. (5) 
J o 



(For equivalence of (|5|) with see \$% .) Here V is the so called Argyres-Kelley projector 

nil 

V... = p B Tr B (...), (6) 
C(t) ... = exp{iC -t}-[H,.. .} exp{-i£ • t}, (7) 

and 
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H = J{a\a 2 + a 2 a{) 



+ 



hu K [g K (ai + a[)(bi K + b\ K ) + G K (a 2 + a\){b 2K + b\ K ) 



(8) 



• Assume, for simplicity, initially factorizable density matrix of the exciton and bath and 
identify p B with the initial density matrix of the bath. 

• Perform the above scaling (including transition to new time t' - we shall, however, continue 
writing t) with also J oc g 2 , 

• Return back from the above 'interaction' picture to the Schrodinger picture. 



(Details of this procedure may be found, for similar models, elsewhere - see, e.g., [41 for a method 
fully following the Davies approach [^|4^j.) Then, after some straightforward algebra, the required 
set of equations for the density matrix of the exciton only (designated as p(t)), that is exact in the 
sense of the above scaling procedure, reads 



d 



( Poo(t) \ 
Pu{t) 

Pia(t) 
\P2i{t) ) 



A B 
C V 



( Poo(i) \ 
Pu(t) 

P22(t) 

Pu{t) 

\p2l(t) J 



Here the blocks 



(9) 



A 



-7T - r T 71 
7T -71 

r T o 




, B 





iJ/h —iJ/h 
—iJ/h iJ/h 



iJ/h -iJ/h 
-iJ/h iJ/h 



V 



i(7| + r i ) + i(e 2 -e 1 )/n 







4(7i + r i )-i(e 2 -e 1 )/n 



(10) 



From the whole set of 9 equations for elements Pij(t), i, j = 0, 1 or 2, we have in @ omitted 
those ones that are separated (the whole set factorizes) and are not important below. Index 
corresponds to the unexcited state where there is no exciton in the system. As for, e.g., the zero 
element in the left lower corner of the square matrix, it is owing (and corresponds) to the above 
stressed importance of the existence of two baths. The following notation has been used: 



7T 



1 



~riTT ^1(9^ k) 2 7771 ~ 



-5{fiuj h 
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n = x^X> K ^) 2 [i + ^X.) - i ]6(huJfi " ei) ' 

Implicitly, we assume everywhere the thermodynamic limit of the bath(s) to be already performed. 
As for Pi and P2, notice that we have assumed the bath to consist of two (sub-)baths designated 
as 1 and 2, each of them connected (forming inherent part thereof) with its own subsystem I or II. 
In order to get rid of the initial condition, i.e. inhomogeneous, terms (otherwise resulting in the 
set ©), we have assumed that the bath is initially statistically independent of the system. We 
impose further condition here: As the two sub-baths do not directly interact, we ascribe both of 
them initial canonical distributions with presumably different initial temperatures T\ and T 2 ■ Then 
Pj = l/(kBTj), j = 1, 2. Unless the opposite gets mentioned explicitly below, however, we shall 
assume the two initial temperatures of the baths equal, i.e. Pi = Pi = P = 1/(/cbT). Then only 
one type of the Bose-Einstein distribution ns(z) = [exp(/3z) — for phonons enters the above 
formulae. Worth mentioning here is also the fact that all the 7's and T's in (|ll|) result as standard 
Golden Rule transition rates (fulfilling the usual detailed balance conditions) between states of a 
localized and absent exciton. No J appears in these formulae. Also this is one of consequences of 
our above scaling procedure that leads to description of relaxation not in the weak-coupling but 
rather intermediate or strong coupling regimes (comparable dephasing which is due to the coupling 
to the bath, and in-phasing which is owing to the above J-term (hopping term)). 

III. RELATION TO THE WEAK-COUPLING LIMIT DYNAMICS 

The weak-coupling limit dynamics is best described in terms of eigenstates of the Hamiltonian 
of the system split off the bath 

2 

H s = ^ e i a \ a 3 + J( a \ a 2 + a\ai). (12) 
i=i 

In our case, one such an eigenstate is known, it is the ground state |0) of the electronic system (no 
exciton in the system). As for the next two states, they read for ei 7^ €2 

\+) = x \l)-<f>\2), |-) = 0|l)+ x |2) (13) 

where 

2Jsign(e 2 - ei) 



2v/(e2-ei) 2 + 4J^[|6 2 - e t \ + y/fe - erf + A.P] 
'2y/(e 2 -e 1 )^ + ^[\e 2 - ei| + ^ - ei )* + 4J 2 ] 



S 



The corresponding eigenenergies (unperturbed by the coupling to the bath) read 



E ± = \ [ £ i + £ 2 T sign(e 2 



ei)V(e2-ei) 2 + 4J2]. 



(15) 



Still for the above case, |+) — ► |1) and |— ) — > |2) when J — > 0. For simplicity, let us consider just 
the case of ei 7^ In the basis of states |0), |+), and |— ), the rigorous weak-coupling dynamics 
1 37, 38] as obtained by scaling (Q) of the exciton density matrix with the above initial condition 
reads 



d_ 

n 



( Poo(i) \ 

P-it) 
Po+(t) 
P+o(t) 
Po-(t) 
P-o(t) 

p+-(t) 
\p- + (t)/ 



( Poo(t) \ 

p++(t) 

p-it) 

Po+(t) 
P+o(t) 
Po-(t) 
P-o(t) 

p+-V) 

\p-+(t)J 



+ 



-0 — T-o 


7o+ 


r ~ 
























f Poo(t) \ 


7+o 


-7o+ 



























P++(t) 


r_ 





-r - 
























P-(t) 











k 


p 


















Po+(t) 
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k* 


















P+o(t) 
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q 












Po-(t) 

















q 


1* 












P-o(t) 























m 









p+-(t) 


























m* 


) 




\p-+(t)J 



Here 



2vr 1 



7+0 = Yn T,(x 2 9l + <p 2 Gl){hu K fn B {huj K )5{E + - huj R ) 



2vr 1 



7o+ = TN + 4> 2 G 2 K ){hu K ) 2 [l + n B (hu K )]5(E + - hu K ), 



9tt 1 

r -o = YN ^ 2 9l + X 2 Gl)(hLorfn B {hu> K )5{E_ - hu K ), 



(16) 



r °- = T ^T,^ 2 9l + x 2 Gl)(hu K ) 2 [i + n B (nu K ))S(E^-nu K ) 



(17) 
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are the Golden Rule transfer rates between eigenstates of the unperturbed Hamiltonian of the 

/ 



system Hs; clearly, e.g. 7+0 — ► 7f when J — > etc. Further, 



EM - VXT ... is the 9 x 9 



V 



diagonal matrix with diagonal elements 0, 0, 0, iE+, — IE+, iE-, —IE-, -%{E- — E+), and 



% : (E+-E-). Finally, 



k = -0.5(7+0 + 70+ + r_ ), 1 = -o.5(r_ + r _ + 7h 



0; 



m = -0.5(70+ +r _), p = 0.5(7+0 + 70+), g = o.5(r - + r_ ). (is) 

Topologically, the weak-coupling relaxation matrix (the square matrix in the second term on the 
right hand side of di~6|)) resembles that in our intermediate or strong-coupling case dH|^). These 
relaxation matrices are written down in different bases, however, and in these bases, they have a 
simple, standard and easily understandable form. Difference between the 'extended' basis |0), |+), 
and |— ) and the 'localized' one |0), |1), and |2) is why the free-motion terms (the first terms on the 
right hand side of ( |l6| ) and (||)) formally differ. Coincidence of the free-motion terms is of course 
complete once they are brought to the same (either localized or extended) basis. This is, on the 
other hand, unlike the relaxation matrix. If we set J — > in all the coefficients in the relaxation 
matrix in the extended basis in the second term in (|l6l), we do not reproduce the relaxation matrix 
of (^) in this basis. Instead, we get, in this way, the relaxation matrix of (|9|) as it is written down 
in (0), i.e. in the localized basis. The point is that 

• in the scaling (^) underlying the weak-coupling regime, the in-phasing (which is due to J 
that is kept constant during the scaling) is automatically, as a consequence of (S) , presumed 
dominating over dephasing processes. Nothing is changed on this feature even if we addition- 
ally take J arbitrarily small. That is why the free-motion term is diagonal in the extended 
basis and the relaxation term describes relaxation in the same basis, i.e. that of eigenstates 
of Hs with J 7^ 0. This relaxation including dephasing is, owing to the form of @, to be 
understood as (in the limiting sense) infinitely slow as compared to the in-phasing processes. 
On the other hand, 

• with our scaling 

g 2 

g — > 0, J — > 0, r — > +00, g 2 T = const, — = const, (19) 

also J is scaled. Hence, as we are allowed, in such a scaling, to keep just second order (in g) 
processes, we must set J = inside all the relaxation superoperator (correction to VC{t)Vp{t) 
on the right hand side of (0)) as the latter is already oc g 2 owing to its proportionality to the 
second power of the coupling to the bath). That is why the relaxation term in (||) describes 
relaxation to eigenstates of Hs\j=o, i.e. in the localized basis. J must be, however, kept 
nonzero in the free-motion term. That is why we get a proper competition between free- 
motion (no transitions in the extended basis, i.e. between eigenstates of Hs\j^o) reflecting 
in-phasing owing to term oc J in <^j, and relaxation going between eigenstates of Hs\j=o, 
i.e. in the localized basis. 
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The situation in our intermediate or strong coupling case connected with scaling (19) thus reminds 
a bit of the Hubbard model when the band and site-local interaction terms are diagonal just in 
the extended and and local bases, respectively. In our case, full equivalence between ( |l6| ) and @ 
appears only in the extreme case of J taken as zero from the very beginning, in both the free-motion 
and relaxation terms. This is because then the extended and localized bases coincide. 



IV. PARTICLE AND ENERGY FLOWS 

Let us return to our above scheme of the slow exciton dynamics as described by (|9HlO|) and let 
us first discuss solution for the case of zero J (for this particular case, correspondence with the 
weak-coupling approach as above can be found quite easily). It is easy to verify that the stationary 
(and in this particular case definitely also equilibrium) solution to (|9|) reads 

Poo = Poo(* -> +00) = 1 - Pn ~ P22, 

exp(-/5ei) exp(-/3e 2 ) , . / r \ <™\ 

P11 = 7= , P22 = 7? , Z = l + exp(-/3ei) + exp(-/?e 2 ). (20) 

As for stationary values of site off-diagonal elements of p, they turn to zero. So as to derive 
(|20|), we have used the detailed balance relations 7f/7j, = exp(— (3e\) and T^/T^ = exp(— ftez)- 
Formulae (|20| ) are still in full agreement with the equilibrium statistical mechanics, in particular 
the canonical distribution. In order to see that, let us realize that we have omitted, for purely 
technical reasons, the two-exciton state with both levels (1 and 2) occupied by excitons. This 
means errors of the order oc exp(— (3{e\ + 62)). Within this accuracy, we can well approximate, e.g., 
the stationary value pn as pn ~ i+cxp7-/?ei) wna ^ * s ^ ne canonical equilibrium probability P\ of 
finding the site 1 occupied by the exciton, as prescribed by the equilibrium statistical mechanics. 
Refraining, on the other hand, for a while from the above omission of the two-exciton state and 
designating the two-exciton state as state 3, we might reconsider the problem on the more general 
level. This would yield the stationary value pn = exp(— j3e\)/Z' and P33 = exp(— (3(e\ + 62)) /Z' 
where Z' = Z + exp(— (3{e\ + £2)). From that, the probability of finding the exciton at site 1 
irrespective of the occupation of site 2 results as P\ = pn + P33 = 1+ex p^ £1 ) Q> m a full 
correspondence with the above reasoning based on discussion of accuracy of our treatment. 

Interesting and important for what follows below is also the dynamics of occupation of site 
1. As we still keep J = 0, we may for this purpose ignore the state 2 at all. Doing so, we shall 
for a while completely neglect 7-f as compared to 7j_. f^] Then the dynamics (time dependence) of 
probability of finding the exciton at site 1 reads 

P 1 (t)=P 1 (t = 0)e-^ t . (21) 



2 This is the standard Fermi-Dirac distribution for excitons that behave as paulions, i.e on-site 
fermions, with zero value of their chemical potential. 

3 It is always 7^/74 = exp(— (3ei) (the detailed balance condition). This ratio is definitely <C 1 for 
ksT -C e\. On the other hand, except in (^T|), this condition is not used below. 
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This corresponds to the probability amplitude of finding the exciton at site 1 in form of e %tt l h "U*/ 2 

whose Fourier transform reads as Lorentzian - T Z\i 2 . , — 7^ . This indicates that the exciton level 

e\ (and similarly for the level €2) is broadened as a Lorentzian with the half- width 7^/2. This is 
what we shall need below. 

All this is very reasonable and known, in different context, for already many years. In what 
follows, we return to a general situation. We shall argue now that the matter will drastically 
change once we put J nonzero. To be more concrete, we are going to argue that under the above 
defined conditions and for still equal initial temperatures of the baths T% = T2 = T, there will be 
a permanent exciton (and also energy) flow between sites 1 and 2. Because of our joining these 
sites with different and mutually non-interacting baths, this implies the existence of the exciton 
(energy) flow also between our subsystems I and II. 

This prediction is perhaps shocking for standardly thinking physicists. In order to deprive 
potential opponents of tempting, at the first sight fully legitimate but still insufficiently justified 
arguments, let us already here state (admit) two important things right here: 

• The weak-coupling theory does not yield the persistent nonzero flows. 

• Application of the canonical distribution to the whole system (consisting of two subsystems, 
each of them having its own bath and electronic, i.e., exciton levels) also yields that these 
flows are in average zero. 

The counter-arguments against such an easy beating off our type of reasoning and results are, 
however, that 

• The weak coupling theory has even potentially no possibility to yield such flows at all. So, it 
cannot serve as an arbitrator. The point is that such flows, as shown below, need a sufficiently 
strong dephasing (as compared with in-phasing processes) to broaden the energy levels e\ 
and 62 (as we shall see below, the transfer is between tails of these levels). The very definition 
of the weak-coupling approach (see (|3|) above) is based on the Van Hove limit leading to a 
negligible role of the dephasing (as compared to the in-phasing), i.e. the approach is unable 
to model the situation we speak about, where we predict the existence (nonzero values) of 
the flows.^ Under the condition of the dominating in-phasing (over the dephasing as in the 
weak-coupling regime), exciton at sites 1 and 2 becomes shared. In other words, a special 
type of a covalent bonding appears between the sites which (as also found in other situations) 
prevents such flows. 

• If we are really right in our prediction that, in the thermodynamic limit, there is a persis- 
tent flow between sites 1 and 2, i.e. between systems I and II, application of the canonical 
distribution is unjustified. The point is that the canonical distribution is based on maxi- 
mizing entropy under solely two constraints: Normalization condition (Trp = 1) and mean 
energy conservation (Tr(Hp) = E =const). If we are right, then at least the third constraint 



4 Notice also that, for, e.g., negligible relativistic corrections and in absence of external magnetic 
field, the canonical density matrix is real in our local basis |1) and |2). Hence, formula (|26| ) below 
yields zero flow between sites 1 and 2 in the canonical equilibrium, irrespective of the Golden-Rule- 
type prediction (|3l|). 
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Tr(Ip) =const (/ being the flow operator) should be added what makes the usual canonical 
distribution improper. Also other approaches used to justify the canonical distribution al- 
ways use, explicitly or implicitly, the ad hoc assumption of non-existence of other persistent 
quantities than energy. So, if we are really right in our arguments here, then, definitely, nor 
the canonical distribution can be used as an arbiter. 

Let us start our reasoning here by deriving a sufficiently reliable form of the exciton flow 
formula. From our Hamiltonian (|l|) and the Liouville equation, we get that 

- — (4ai) = ^{[a{a 1} H]) 
= -^-3?m(4ai) + 9K^K^rn( b XK a i)- ( 22 ) 



The last mean value can be calculated, using principle of the adiabatic switching on the interactions, 
as 

W\ K a 1 ) = j t (b[ K a 1 ) = ~([b\ K a 1 ,H}), <5 ^ + . (23) 

Here, terms containing J should already be omitted if we work to the second order in our small 
parameter g only (remember that J oc g 2 - see (|l9|)). Thus, within this accuracy, 

(b{ K ax) -=g K hw K {n B (tiu>K)[l - (a! a i) ~ t 1 + n B (huj K )}(a\ai)}. (24) 

n(jj K — e± + iho y/Jy 

So, fl22] ) reads, within the required accuracy, as 

As the last two terms express the exciton number imbalancing owing to transfers 1 <-» 0, the proper 
formula for the real 1 <-> 2 flow is connected with the first term on the right hand side of (|25|). 
Thus, the 1 <-> 2 exciton flow (taken as positive if flowing from 1 to 2) reads 

/ = ^-3?mpi2. (26) 

The fact that I is determined by the (imaginary part of the) site off-diagonal elements of the 
particle density matrix follows already from the elementary quantum mechanics. 

The long-time (stationary) value of the p\2 element of the density matrix can be found, however, 
from (||) (by setting the time-derivatives zero), incorporating also the normalization condition 

3=0 

After a simple algebra, the result is 



/= 2vr 2 1 |(7j+£i) 
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7 i r i + 7i r T + 7T r i + x[ 7i + r i + 2 7T + 2r T ] ^ v ; 

Here 

X = ^j2 1 §(7i+ri) 

Since we are obliged to stick to the required accuracy, we can deal, using the formalism allowed, 
with just the leading terms. Hence, we shall omit the terms oc X in the denominator assuming 
that 

X < llT l ■ (30) 
7|+T| 

This means that, up to terms of higher than sixth order in g (still remember that J oc g 2 ), 



7t r l ~ 7|F T 

7^ + 7i r T + 7T r ; 



_2tt 2 1 ICTi + Ij) 

~ X J '^7^7 , -p \io i r 7^[Pn-p22\- (31) 

71 7r [f(7j. + rj.)] 2 + [e2-ei] 2 

Interesting point is that the expression for the exciton flow on the right hand side of ([H]) is correct 
even without assuming (|30|). This follows from (|28|) by taking into account that from (^), we obtain 
the asymptotic-time populations 

= 7T r i + x( 7T + r T ) 

Pn 7i r i + 7i r T + 7T r 4 + x( 7i + r i + 2 7T + 2r T ) ~ z' 

P22 = 7i r T + x (7T + r T ) 

7ir i + 7i r T + 7T r i + x( 7i + r i + 2 7T + 2r T ) z v ; 



In the last approximate expressions, we have again used condition (|3fj|). 

Clearly, expression (|3l]) for the exciton flow is, quite surprisingly at the first sight, clearly 
nonzero. Already this is remarkable as we have to realize again that the exciton transfers energy 
and the transfer channel l<->2 is the only channel connecting our systems I and II and able, within 
our model, to transfer energy between them. Before getting into more physical details connected 
with this observation, let us also comment that expression (|3~l"D 



has a proper total balance structure with transitions 1 — > 2 and 2 — ► 1 (contributing to (31) 
by terms oc p\\ and oc p22> respectively), and 

is fully compatible with (in fact, it is exactly equal to) the second-order Golden Rule of 
quantum mechanics (J is the matrix element of the transfer part of the Hamiltonian between 
states of the exciton at sites 1 and 2) assuming that the energy conservation law is properly 
broadened owing to exciton decay processes. 
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This broadening means that the exciton transfer is neither at level e\ nor at level £2(7^ ei) but 
generally at arbitrary energy in tails of the two broadened exciton levels. (Realize that the exciton 
is, as generally in nature and as also anticipated in our model, just a finite life-time quasiparticle.) 
This interpretation is clearly confirmed by the fact that one can rewrite (j31[) also as 



2* 2 f w 1 'ill 1 |£i 

and the fact that any quasiparticle exponentially damped with the decay rate 7, i.e. the survival 
probability amplitude 

(ai(i)4) = exp(-iei • t/h - -yt/2), (34) 

has its energy level (here ei) broadened into a Lorentzian with the energy half-with hj/2. See also 
a comment in this respect above. The forms of ([H|) and ( |33"D thus leave only very limited space 
for speculations about validity of our approach. 

Our statement about persistent energy transfer between subsystems I and II kept at (initially) 
equal temperatures is already justified by our expressions for the exciton flow /. Nevertheless, let 
us raise the (in a way) subsidiary question what is the energy flow between the two subsystems. 
According to the above quasiparticle interpretation, one would expect that the energy flow is 

n 27r 72r i f + °° 1 hi 1 f r i A 

Q = ~r J [Pn- P22 - - 71 19 , f r- ■ - h -ede 

h J-00 vr [§ 7i ]2 + [ e _ £l ]2 vr [| r j2 + [ e _ £2 ]2 

= -r-J [pn ~ P22] ■ —Tn ( , r x 12 , r 12 ' ^ 

The problem is with general justification of this formula. In fact, one should define the energy flow 
Q between systems I and II in full generality as 



or equivalently 



Q=-j t m (36) 



Q = J t (Hli). (37) 



Here Hi and Hjj should have the meaning of energies of the subsystems I and II, such, that 
Hj + Hj[ = H (in order to have ( p6| ) compatible with (p7|)). Because J must be assumed nonzero, 
these evidently cannot be the Hamiltonians Hj and Hjj introduced in (|l|). Really, making this or 
any other trivial identification of Hj and Hn leads to hardly interpretable results. The physical 
reason for that is that nonzero values of J cause effects like J-dependent renormalization of the 
exciton coupling to the bath. Its exact form is unknown so that one can ignore it only when the 
corresponding coupling constants are negligibly small (when there is nothing to be renormalized). 
That is why, for very small g K (i.e. negligible 7 j_), one can define the energy flow properly and 
reliably by (|36|), identifying (in this particular case) Hj with Hf. Then ( |36| ) reduces to the above 
formula (|35| ) with negligible 7 |. In the opposite limiting case, when G K gets very small (i.e. with 
negligible T^), ([37]) also reduces to the above formula (|35[) , this time with negligible 1^, provided 
we identify Hjj with Hjj. 
So summarizing: 
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• Suggested interpolation formula for the energy flow (|35|) can be properly justified in the two 
above mentioned limiting cases. Then it definitely yields nonzero values of the energy flow 
between our subsystems when e\ ^ ei and still T\ = T2 (equal initial temperatures of the 
two baths). 

• Formula for the exciton flow (|3l[), or its equivalent form (p3|), can be, in the above way, 
properly justified for all the values of the parameters involved. It always, for e\ 7^ €2 and 
T\ = T2, yields nonzero exciton flow between the subsystems, always in the direction from 
the system with higher population to that one with the lower population of the corresponding 
exciton level. As the excitons bear energy, this clearly implies (and fully corresponds to the 
above) nonzero mean energy flow. 

The last point to be stressed here results from rather a trivial reasoning based on the above 
formulae and continuity, with respect to, e.g., temperature T2, of the results obtained. On the 
other hand, it is of highest importance from the point of view of interpretation: It is easy to see 
that if we take (as always here) energies of the exciton levels in our two systems different, i.e. 
e i 7^ £ 2j but allow the initial temperatures of the two subsystems be different {T\ 7^ T2), we can, 
for any values of e\ 7^ €2 and arbitrary T\, find such a value of the initial temperature T2 of the 
second subsystems that 'separate-equilibrium' (for J = 0, i.e. for separated subsystems I and II) 
values of population of the disconnected exciton levels e\ and €2 become equal. In other words, 
P11 = P22- According to what has been argued above, however, this makes, upon reintroducing the 
(sufficiently weak - see (|30|)) exciton transfer channel (setting J nonzero though small enough in 
the sense of (|30"|)) both the exciton / and the energy flows Q (simultaneously) zero - see formulae 
(HD, flU, and (||). 

In particular, assume now that, e.g., 62 > e\. Then clearly, arguing just from the above 
'separate-equilibrium' values of p\\ and p22> the temperature T2 for which the values of p\\ and 
P22 get equal, is greater that T\. Let us call this temperature T^ lt - It is easy to show that 
rpcnt _ T\t2/ e l- Clearly, for temperatures T2 greater than T\ but less than , the energy 
(exciton) flow is nonzero and still going in the direction from the subsystem I with the higher 
exciton population (lower temperature T\ < T2 but essentially less excitation energy ei < 62) to 
that one with the lower exciton concentration - subsystem II. This means, the exciton and energy 
flow goes, fully surprisingly, against the temperature step. This conclusion is of highest importance. 



V. TOWARDS THE SECOND LAW OF THERMODYNAMICS 

There are several formulations of the second law. The form by Clausius from 1865 involves 
entropy that was not discussed here. Existence of entropy is, however, in fact consequence of 
three main formulations of the second law, that one by Thomson (1849), Clausius (1850), and 
Caratheodory (1909) (for connections to entropy see || or standard textbooks on axiomatic ther- 
modynamics). The statements are (cited according to ||): 

Thomson (as Baron Kelvin of Largs since 1892) |42(| : No process is possible, the sole result of 
which is that a body is cooled and work is done. 

Clausius [|43||: No process is possible the sole result of which is that the heat is transferred from, 
a body to a hotter one. 
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Caratheodory [44|: In any neighbourhood of any state there are states that cannot be reached 
from it by an adiabatic process. 



The words 'sole' imply in particular that 

• in the Thomson formulation, the process should be cyclic, without any compensation (ad- 
ditional heat transfer to another and cooler body). A (thought) machine working in such a 
style is often called 'perpetuum mobile of the second kind'; 

• in the Clausius formulation, the process is not aided from outside. 

Significance of this law and consequences of its potential violation were perhaps best described 
in the Introduction of §. For 150 years, nobody really questioned this statement based on un- 
countable number of observations from our everyday life. One should add and stress therefore that 
our thought (i.e. still not real in the sense of in-nature-existing) systems governed by quantum 
mechanics and comprising macroscopic baths are macroscopic, i.e. they should obey the Clau- 
sius (and not only this) formulation of the second law provided that the quantum mechanics and 
thermodynamics are always, including the quantum- and at least the macro- world, compatible. 
Sometimes, even absolutistic statements in favour of the unconditional validity the second law in 
the macroworld stemming from our everyday experience appear ("...No exception to the second 
law of thermodynamics has ever been found - not even a tiny one..." [jig]). As our conclusion at 
the end of the previous section shows, the opposite is true in our model here. One should mention 



here the 'pawl and ratchet' systems originally suggested by Feynman [23 which are also often cited 
also in connection with the second law. These systems, however, so far fail in practical attempts 
to violate the second law [Q, in full agreement with the Feynman |23| analysis. On the contrary, 
behaviour of the above model does, as argued above, contradict the second law. It is of course not 
for the first time that such a mathematically well-founded behaviour incompatible with standard 
thermodynamics appears - see [TT|, 8| and, for other models, papers cited above or therein. The 
present model is, on the other hand, perhaps the simplest one. Worth mentioning is also that the 
above criteria for the energy flow against the temperature step may easily be compatible with, e.g., 
even room or higher temperatures. This is unlike, e.g., [14]. In the infinite temperature (i.e. in the 
classical) limit, however, the effect disappears. 



VI. TOWARDS THE ZEROTH LAW OF THERMODYNAMICS 

In order to be specific, let us state what this law (so often, especially in the mechanical context, 
understood as trivial) says: // system A is in equilibrium with systems B and C then B is in 
equilibrium with C (see, e.g., |}4q| ). It helps to introduce thermodynamic temperature, chemical 
potential etc. Though universal validity of this law has already been questioned (as far as its 
form for equality of chemical potential of one sort of species in different phases in equilibrium is 
concerned) - see |2£| or in the implicit form in 0], no special attention has so far been paid to this 
fact. That is why we should address the question, in connection with our model above, again. 

Let us fix, in the above model, again the situation with ei < ei and be T\ arbitrary positive. 
Then clearly T^" 4 = T^/ei > T\ and we set the initial temperature of the second bath T2 = T§ . 
Hence, upon establishing a contact between subsystems I and II by taking J slightly (in the sense 
of ([30])) nonzero, we get no energy or exciton flow. Now, we can invoke standard thermodynamic 
definition of what it means to say that two bodies in a contact are in mutual equilibrium. The 
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definition reads that introducing arbitrary obstacles hindering flows between the bodies does not 
change their thermodynamic state. In this sense, this is exactly the situation we are now in: We 
have two bodies in a contact where there are no flows between them. Hence, the thermodynamic 
state cannot be violated by any obstacle setting these flows zero and not influencing otherwise the 
state of the systems because the flows are already zero. (Notice that no other flows but those of 
exciton or energy can exist in model.) Remind, on the other hand, that we have such a strange 
thermodynamic equilibrium (in the above thermodynamic sense) that temperatures of both the 
systems (those of their baths) are different. This can clearly lead to other contradictions with the 
standard thermodynamics as we are now going to show. 

Assume now that we have still another (third) exciton level and still another (third) thermal 
bath attached to it. In other words, we complement our Hamiltonian by terms 

AH = e 3 a\a 3 + V huj K b\ K b 3K + — L ^ h K huj K (a 3 + 4)(6 3k + &|j 

ViV _ 



+ ~7tt H K hoj K (ala 3 + 4a 2 )(&2K + h K + b\ K + b\ K ) + K{a\a x + a\a 3 ). (38) 



Clearly, the last two terms on the right hand side of (^) provide interaction of subsystems II and 
III (this type of the phonon-assisted interaction exists just in the diabatic (non-rigid) basis |fl7f ) 
and that of the subsystems III and I (the latter interaction and the induced exciton transfer is 
for simplicity assumed coherent, like that one of the subsystems I and II). Assume also that the 
exciton energy £3 equals to that of the exciton in subsystem I, i.e. £3 = e\. Next, we assume 
that again, the initial density matrix is factorizable into a product of density matrices of all the 
subsystems, so that there are no exciton-bath initial statistical correlations between any two of the 
three subsystems. Finally, assume that the density matrices of all the baths are initially canonical. 
Corresponding temperatures are assumed as T 3 = T 2 > T\ (the last inequality being assumed 
already above). 

Let us for a while set J = K = 0. Then we have our subsystem I fully separated and the 
dynamics goes between subsystems II and III only. Exactly in the same way as above (i.e. using 
the same type of scaling), we get a closed set of equations for the matrix elements of the exciton 
system. This time, however, the situation is simpler as compared to that above. First, our coupling 
between subsystems II and III is assumed as bath-assisted. In connection with that, the set of 
equations for the site diagonal as well as site off-diagonal matrix elements of the exciton density 
matrix factorizes so that equations comprising the diagonal elements contain only the diagonal 
elements (that get separated from the set for the off-diagonal elements), reducing in form to the 
Pauli master equations |48|| . Properties of these equations are sufficiently known. So we shall not 
repeat the calculations and refer the interested reader to any elementary textbook of kinetic theory. 
The result for the asymptotic exciton occupation probabilities reads as in the standard equilibrium 
statistical thermodynamics, i.e. 

exp(-/3 2 e 2 ) 



P22 



1 + exp(-/3 2 e 2 ) + exp(-/3 2 ei) ' 



exp(-/3 2 ei) „ 1 , . 

1 + exp(-/^ 2 e 2 ) + exp(-/? 2 ei) k B T 2 



Let us stress the following points: 
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• Values ( |39| ) properly reproduce, within our approximation (omission of multi-exciton states, 



i.e. with errors oc exp(— fofa + standard equilibrium statistical mean exciton numbers 
at sites 2 and 3, i.e. l/[exp(/32£2) + 1] and l/[exp(/?2ei) + 1] (see a comment above concerning 
appearance of these Fermi-Dirac distributions for excitons). In fact, as already argued above, 
reintroducing the multiple-exciton states would reproduce these values exactly. These Fermi- 
Dirac distributions are, on the other hand, proper mean number of excitons at sites 2 and 3 
for H K = 0, i.e. separated subsystems II and III. Hence, establishing or cancelling the above 
contact between the latter two subsystems does not change the stationary (equilibrium) 
exciton populations at the corresponding sites. The same may be shown to apply to phonon 
populations in the corresponding baths. 

• Assume now H K ^ 0. Owing to the incoherent (bath-assisted) character of the above coupling 
between subsystems II and III, different populations of levels 2 and 3 (in accordance with 
standard statistical thermodynamics) do not contradict the fact that (exciton mediated) flows 
between subsystems II and III remain in equilibrium exactly zero. This is due to the fact 
that real 2 — > 3 and 3 — > 2 transition rates are proportional to P22 x {1 + l/[exp(/?2ftk> K ) — 1]} 
and x l/[exp( / 92^w K ) — 1]}, respectively. Here, by the energy conservation law, hjj K = 
62 — ei > 0. The multiplicative factors at P22 and P33 are phonon statistical factors describing 
phonon-assisted induced as well as spontaneous processes. So, the transfer rates 2^3 and 
3 — > 2 are in equilibrium exactly equal, mutually cancelling their contribution to the exciton 
as well as energy flow between subsystems II and III (transferred energy is ei + fko K = £2). 
Here, for simplicity, we have assumed (in accordance with assumptions underlying validity 
of the Pauli equations) that |e2 — £3!/^ = |£2 — £i|/^ is appreciably greater than the sum 
of broadenings of levels 2 and 3, i.e. that the transitions are practically (exciton+phonon) 
energy conserving. These arguments are what underlies the detailed balance conditions in 
the Pauli master equation theories yielding the same conclusion. 

These are the characteristics of the mutual equilibrium (according to the above thermodynamic 
definition) state of subsystems II and III, and also of the internal equilibrium states of the isolated 
subsystems II and III taken separately. With this states, let us now put H K = (we split subsystems 
II and III), keep J = but put K ^ 0. Let us repeat: We have the two subsystems (I and III) 
with equal exciton energies e± = £3 (this case may be treated as a limit e\ — £3 — > in the above 
formulae), initially in internally canonical states of both the subsystems (and their baths), with the 
respective temperatures T\ < T3. In accordance with what has been said above about development 
of two such subsystems (subsystems I and II above) with their respective baths, the asymptotic 
(stationary) populations of the exciton levels 1 and 3 only slightly change upon establishing contact 
between the subsystems. Asymptotically, they read 




exp(-/3i£i) 




exp(-/? 2 £i) 



1 



(40) 



Clearly, because fa < the populations pn and P22 are different. So, according to (pl|), (|33|) as 
well as (|35|), there are flows between the subsystems I and III, i.e. we have no equilibrium in the 
thermodynamic sense. 
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Thus, summarizing, we have subsystems I, II and III which all have well defined temperatures. 
Upon establishing just the above specific contact between I and II, the systems stay, in the sense of 
the thermodynamic definition, in equilibrium. Similarly, upon establishing just the above contact 
between subsystems II and III, the thermodynamic equilibrium is not violated. Thus, the zeroth 
law of thermodynamics should be, according to thermodynamics, well applicable. It states that 
establishing contact between subsystems I and III should preserve their mutual equilibrium state. 
As seen above, however, the opposite is true. 

VII. CONCLUSIONS 

We have investigated one standard and rather trivial model that allows rigorous treatment 
by methods of the quantum theory of open systems. The obtained behaviour contradicts what is 
prescribed by the second as well as zeroth laws of thermodynamics. In connection with previously 
expressed doubts about universal validity of the second law in specific situations, this extends 
challenges to general compatibility of such two basic scientific disciplines as the thermodynamics 
and the quantum theory. 
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